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1 Introduction 


The problem of synthesis with bounded controls (SBC) of the time opti¬ 
mal control [T] is one of the well-known optimal control problems. The 
problem of SBC is stated as follows: 

Given a system of differential equations, 


x = f{x,u), X € M", u £ C M’’, 0 £ int fl. (1.1) 

It is required to find a control u of the form u = u{x), with values in the 
set such that the trajectory of the system 


X = f{x,u{x)), 


( 1 . 2 ) 


starting at an arbitrary point Xq, terminates at a given point Xi in min¬ 
imal time. 

In this case, the control u{x) is time optimal and satisfies the Bellman 
equation [2]: 

") = -'■ (i-s) 

where T(x) is the time of motion along the trajectory of system (II. 2|) . 
which corresponds to the control u{x). Let us denote by r(x)| fYj| , 
r(x)|jY^ the time derivative along system (II.1|) . (11.21) respectively. 
Hence equality (11.31) takes the form 

min r(x)|{X^ = r(x)|{X2J = “1’ 


which means that the derivative by virtue of system dOJ) of the time 
optimal T{x) from an arbitrary point x to a given point xi is equal to 
— 1. Obviously this equality holds at points where the derivative exists. 

If we renounce the time optimality, we would consider the admissible 
positional synthesis problem, which consists of the construction of the 
positional control u = u{x), which satisfies a given restriction, i.e. u £ H. 
As a result, the trajectory of system starting at an arbitrary point 

Xq, terminates at a given point xi at some finite time T{xo). Furthermore, 
we assume that xi = 0 and /(0,0) = 0. 

For solving this problem, the method of the controllability function 
was introduced by V.I. Korobov [3]. This method is based in the con¬ 
struction of the controllability function 0(x) (0(x) > 0 for x 0 and 
0(0) = 0), as well as, positional control u{x) = u(x,0(x)), such that 
inequality 


E 


dO(x) 

dxi 


fi{x,u{x)) <-/3Q^ ^/“(x) 


(1.4) 
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is satisfied for some /3 > 0, a > 0. The realization of this condition 
guarantees that the trajectory arrives to the origin at finite time. We as¬ 
sume that the controllability function 0(a;) is continuously differentiable 
at a; A 0. In case that inequality (11.41) holds for a = oo, the function 
0(a;) is a Lyapunov function. 

In [4] a general method of constructing a controllability function and 
a control, which is solution of the admissible positional synthesis, was 
introduced for the linear control system 

X = Ax + Bu, X G M", U&W, (1.5) 

where A and B are (nxn) and (nxr) constant matrices, respectively. In 
particular, the function (p{s) = 1—s for s G [0,1], (p{s) = 0 for s > 1 
generates a controllability function Qtp{x), which is the time of motion 
IE3 with the control u^{x) from an arbitrary point x to the origin 
according to the system dia. 

We are interested in finding a wider set of pairs of functions: the 
function of controllability 0(x), which is the time of motion and the 
control which solves the synthesis problem. 

In the first part of the work, we consider in detail this problem for 
the canonical system. The second part is devoted to the solution of 
the positional synthesis control problem for the completely controllable 
system with a restriction of the form u G fl = {u : ||u|| < d}, 

d > 0. The controllability functions 0(x), which are the time of motion 
from the initial point x to the origin, satisfy the equality T{x) = 0(x), 
and the derivative of the controllability function by virtue of the system 

X = Ax -I- Bu{x) (1.6) 

is equal to —1: 0(x)| JY^ = —I. 

In the present work, we use the controllability function method which 
was first introduced in [3]. 

2 Preliminary results 

We consider the canonical control system 

x = Aox + bou, ( 2 . 1 ) 

where 



/ 

0 

0 . 

. 0 

0 

0 



/ 

1 

\ 

Ao = 


1 

0 . 

. 0 

0 

0 

, bo = 


0 




1 

0 

0 . 

. 0 

1 

0 



1 

0 

/ 
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with restrictions on control |m| < d. We will choose a control as a function 
depending on states coordinates u = u(x) such that the trivial solution of 
system for u = u{x) would be asymptotically stable. We will refer 
to this control as the auxiliary control. Real numbers oi, 02 ,..., On will 
be chosen in such a way that A" — — ... — oi = 0 would have roots 

n 

with the negative real part. In addition, set u(x) = ^ atXi = {a,x), 

i—1 

where a = (ai,..., a„)*. In this case, system (12.1|) has the form 

X = Aix, (2.2) 

where 

^ Oi 02 ... ari-2 O-n-l On \ 

, 1 0 ... 0 00 

All = ; 

\ 0 0 ... 0 I 0 y 

its trivial solution is asymptotically stable. In this way, the auxiliary 
control solves the stabilizability problem for system m in all the phase 
space, but does not satisfy the given restriction. Since Ai is stable, 
we can find a positive definite quadratic form V{x) = {Fx,x) which 
is a Lyapunov function. The derivative of V by virtue of system (123) 
represents a negative definite quadratic form (—Wx,x), with positive 
definite matrices F and W. Since 

— {Fx, x) = {Fx, x) + {Fx, x) = ((T'Ai + A*F)x, x), 
dt 

this equation is reduced then to the Lyapunov matrix equation 

FAi + AlF = -W, 

which for given matrices Ai and W has a unique solution. 

Let be a positive definite matrix, the form of which will be specified 
later. Denote by D(0), H diagonal matrices of the form 

/ 2i-l \n / 2? — 1 \ ” 

D(0)=diag(0- —iJ = diag(^- —j . 

Let oq be a positive number which will be determined later. For 
each X 0, define the controllability function 0(x) as a solution of the 
following equation: 


2ao0 = (F>(0)FF>(0)x, x). (2.3) 
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Set 0(0) = 0. This equality and equation (12.31) determine the function 
0(a;) which is continuous for all x and continuously differentiable for 

X A 0- 

We will look for control u{x) in the form 


n 

u{x) = e-i{x)a*D{e{x))x = 


(2.4) 


where F and the vector a will be selected in a such a way that the 
controllability function 0(x) must be the time of motion from the point 
X to the origin. In other words, the time derivative 0(x) by virtue of 
system (EB with the control u{x) (El satisfies equality 

0(x) = -l. (2.5) 


Lemma 2.1. Let the controllability function 0(x) satisfy equality V2.5]) . 
Thus, 


Oi = - 


n(n + 1) 
2 


( 2 . 6 ) 


Proof Set y{Q, x) = D{Q)x. The controllability function 0(x) for x A 0 
then satisfies equality 

2ao0(x) = (Ay(0(x), x), y(0(x), x)), (2.7) 

and control (12.41) has the form 

u(x) = 0~^ (x)a*y(0(x), x). (2-8) 


In view of equalities 

D(0)AoD-^(0) = 0-Uo, D(0)bo = 0-ho 

the derivative of the function y(0(x),x) by virtue of the system (12.11) 
with this control has the form 

y(0(x),x) = 0~^(x) (^0(x)ff + Ao + boa*J y(0(x), x). 

Therefore, from equality (El we find that the derivative of the control¬ 
lability function by virtue of system El with control u(x) (12.81) is given 
by equality 

• _ ((F(Ao + bpa*) + (Aq boa*)*F)y{0{x),x),y{0{x), x)) 

{iF-HF-FH)y{0{x),x),y{0ix),x)) 


Translation of the paper: Mathematical Physics, Analysis and Geometry, (Russian) 2004, Vol.ll No.2 


212 






Controllability function as time of motion. I 


Hence, by the conditions of the lemma we have 

^0 + b^a* + —I — + bf^a* + —I — F = 0, (2.9) 

where I is the identity matrix. From this equality, we conclude that 
the matrix F^i^Ao + boa* + — H)F~^ is a skew-symmetric matrix, 

therefore its characteristic equation has the form 

det (f^ -I- boa* ~ ~ 




= (deti^-5)2det F Ao + boa* + -I - H ] - XF] = 0 


which is equivalent to 

det ( ( Ho -|- boa* 


I-H] -XI] = 0. 


( 2 . 10 ) 


In addition, it has eigenvalues with zero real part. Let us write equality 
(12.101) in the following form: 


/ 

0 = det 


A — 1 — tti 
-1 


-a2 

A-2 


On—1 
0 


C-n '\ 
0 


\ 0 0 


— 1 X — n j 


n n— 1 n 

= ]^(A-j)-^aj (A-i)-a„ = A"-(ai-f n(n-k l)/2) A”“^-.... 

1=1 i=i i=i-l-i 


From the fact that the roots of the last equality have zero real part and 
are complex conjugate, we have that Oi = —n(n -I-1)/2. 


Lemma 2.2. Let F = be a positive definite matrix. 

Thus, the matrix 



( i 

fl2 ■ 

■ fin 

\ 

G{0 = 

/21 

/22 ■ 

■ f2n 



V fnl 

fn2 

fnn 

/ 


is a positive definite matrix and det G(^o) = 0 for f > fo- 
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Proof. The determinant det G(^) is a linear function on f, i.e., 
G(C) = + b, hence by Sylvester criteria a > 0. Therefore, if 

det G(/ii) = a/ll + 6 > 0, then for any ^ > /n we obtain that det G(^). 
Moreover, there exists a unique such that det G(^o) = O'fo + ^ = 0, 
and obviously, for auy ^ we have that det G{^) > 0. 


Lemma 2.3. 

i) The (sxs) matrices 

Ps.k=( I ) 

are positive definite ones. 

ii) For every fc £ N, we have 



fc+i 


fe|i 


fe +2 


\ s+/c—1 s-\-k 


s+fc-l \ 
s+k 

2s+fc-2 / 


As,fc = det 


( - - - 

\i+ j + k - 2 


1 

i + j + k — 1 


> 0 . 


fc G N, 


Hi) 


The determinant of the (nxn) matrix P of the form 


/1- 

n(n+l) 

1 

1 

\ 

2 

2-3 

’ ■ * n(n+l) 


-1 

1 

1 



3-4 

■ ■ • (n+l)(n+2) 



0 

1 

1 



4-5 

(n+2)(n+3) 


V 

0 

1 

1 

) 

(n+l)(r!,+2) 

(2n-l)2n 


is equal to zero, and its rank is equal to (n—1). 


( 2 . 11 ) 


Proof, i) The positive definiteness of the matrices Pg.fe follows from the 

positive definiteness of the Hilbert matrix ( j 

V ^ / i,j—l 

The proof of this fact readily follows from the integral representation 


P.t.k = 


1 

( ^ \ 

/ 

t 

J 

0 

\ / 


(l,t,...,C )d 


The matrix Pg^k is then positive definite one [5]. 

ii) Let us represent the determinant Ag ^ in the form 


As.fc = det 


i j k — 2 i+j + fc—1 
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Add to the i-th row {i > 1) all the next rows, that is, the rows with 
numbers *+l, ..., s, and extract from every z-nth row (z = 1,..., s) of 
the obtained determinant the value (s—z+1) and from every j-th column 
the value l/(s+j+fc—1) j = 1,..., s. We obtain 




{s + k)... (2s + k — 1) 


detPs,k- 


From there and in view of i), it follows that Ag^k > 0, fc € N. 
iii) Let us write the determinant of the matrix P in the form 


1 - 


i(n+l) 


-1 


n+1 


3 

4 

' ' n+1 

n+2 

1 

1 

1 

1 

(n+1) 

n+2 

2n-l 

2n 


( 2 . 12 ) 


We transform this determinant in two different ways. In the beginning, 
we add to every j-th column (j > 2) of this determinant all the next 
ones; the columns with numbers j+1, ..., n. We have 


detP = 


1 - 


n(n+l) 

2 

-1 

0 


(n + 1)... 2n 


n—1 
2(n+l) 
n—1 
3(n+2} 

n—1 

(n+l)2n 
1) (i_ 

— {n + 2) 
0 


"-i)++i) 

n(n+l) 


n(n+2) 

(n+l)(n+2) 

= 

2 

1 


(2n—2)2n 

(2n-l)2n 


1 

1 1 


2 

n—1 n 


1 

1 1 


3 

n n+1 

1 

1 1 


n+1 

2n-2 2n- 

1 


(2.13) 

On the other side, by adding the rows with numbers from z + 1 till n to 
the z-th row (z > 1) of determinant ()2.12l) . we obtain 


det P = 


n(n+l) 

2 

-1 


(n + 2)... 2n 


2(n+2) 
n—1 
3(n+2) 

1 

{n+l)(n+2) 

_ n+1 

2 ' 
1 

n—1 

0 


(n—l)2n n(2n) 

n—1 n—1 

n(2n) 

1 


(n+l)2n 

_ 1 

(2n-2)(2n-l) (2n-l)2n 

1 1 


n—1 
1 


1 

n+1 


0 


1 

n+1 


2n—2 2n—1 


(2.14) 
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Denote by Bi and B 2 the complementary minors with numbers (1,1) 
and (1,2), respectively, which appear in the RHS of (12.131) and (12.141) . 
Thus 


det P = 1 — 


i{n + 1) 




(n- 1)! 


n + 1 J (n + 1) • • • 2n ’ 


detP=|-d:i±hB, 


B, 


(n-iy. 


— 1 / (n + 1) • • • 2n ’ 


and 


Pi = 


- 1 


Bo, detP = 0. 


Since Pi A 0; then rank(P) = n — 1. The Lemma is proved. 

We determine matrix P, which appears in equation (12.31) . and the 
vector a, which appears in control of form (12.41) . From equality (12.9L we 
have 


Aq + boa* + —I — H 


P- 


P- 


boa* 


-H 


= 0 . 

(2.15) 


Denote P„ = diag ((—l)*“^/(i — 1)!)”_^ ■ We must hnd the matrix 
P~^ with the form P~^ = DnCDn, where the matrix C is a Hankel 
matrix C = (ci+j)"^.^^. We find the matrix C. To this end we write 
equality (I2.15P in the form 


_/ _ H+D-^AoD^ + D-Hoa*Dr 


C+ 


+c{^I-H+D-^ADr, + D-%a*D^ =0. (2.16) 

Denote by Q = [qij )^the matrix in the LHS of equality (I2.16|) and 
write it in the following form 

n — ( \ 

^ V Q 21 Q 22 j ’ 

where Qn = qn, Qyi = (gi 2 , ■ • ■, Qin), <521 = Q 121 Q 22 = (9o)"j=2- 
Hence the elements of the matrix Q are given by the following equality: 




n—1 

<711 = 2(1 + ai)co + 2 Oj-i-iCj, 

i=i 

n—1 

qu = -(* - l)cj-2 + (1 +di + i)cj_i + '^'dj+iCi+j-i, i = 2 ,...,n, 

i=i 

qij = -{i+j - 2)c^+j-3 + {i+j)ci+j-2, i = 2 ,...,n, j = 2 , ... ,n, 

(2.17) 
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where ai, ..., a„ are the components of the vector a = DnU. Since Q 
is a zero matrix, then from the equalities Q 22 = 0 we obtain 


(2n — l)2n 

= (j + l)0- + 2)"^"-^ 
and equality (12.171) has de form 


j = 1,... ,2n - 3, 


(1 + ai)co «j+i ^ „ 

2n(2n - l)c 2 ri -2 ^(j + l)(j + 2) 

n . ~ n-1 ~ 

"2n(2n-l)c2n-2 2-3 ^^(j + 2)(j + 3) 

_ 1 4 + ai Qj+i ^ f. 

3 3-4 +^^0- + 3)(j+4) ’ 


(2.18) 


1 111— — 

1 ^n + l + ai ^ ®j+i _ Q 

n n{n + l) ^ (j + n)(j + n + 1) 

where di = oi is defined by equality (1^ . Consider equality (12.181) . as 
a linear system of equation with the respect to the unknowns 
_Co_ 

(2n — l)2nc2n-2 ’ 

which in vector form can be written as 


Py = yo, 


(2.19) 


where P is a (nxn) matrix defined by equality (I2.11|) . In addition, y, yo 
are vectors of the form 


y= 


Co 


(2n — l)2nc2„-2 


, 0-2 ^ ^ CLj 


yo= 0,- 


Oi 


2-3 


Ql 

'3-4 


Oi 


i{n + 1) 


This system represents a subsystem, since the matrix C = 

and the vector a* = —^hQD~^C~^D~^ satisfy the equality (12.161) . Ac¬ 
tually, the equality (12.161) for such a selection of the vector a has the 
form 

-I — H + D„ C + c(-I — H + D„ ^ = bnbn. 


This equality holds if one inserts the matrix C in the last equality. 

Therefore, rank(P) = rank(P,yo); by virtue of iii) and Lemma 12.31 
we have that rank (P) = n — 1. 
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Let us find all the solutions of system (12.191) . Consider matrix 
((n—1)X(n—1)), which has the form P 


P = 


I -1 

0 


1 

3- 4 
1 

4- 5 


n(n+l) 

{n+l)(n+2) 




n 1 1 / 

\ " (ra+l)(n+2) (2n-2)(2n-l) / 


This matrix, by because of ii) of Lemma 12.31 for k = A and s = n—2, is a 
nondegenerate matrix. Denote by = and! -\- d", where 


d'=[- 


1 


1 


(n + l)(n + 2)’ (n + 2)(n + 3)’ ’ (2n — l)2n 

^„_/_3 + ai oi oi 


2-3 ’ 3-4’ ■■■’ n(n + l)^ 

Furthermore, by y we denote the (n—l)-dimensional vector 


y = 


Co 


(2n - l)2nc2„-2 


, 02, ..., a„_i 


and consider a system of equations Py = yo with respect to y. This 
system has a unique solution y = 


Co 


(2n — l) 2 nc 2 n -2 
1 


= 1(a;5„ + a") = ^(^ 


,(-l) 


n —1 


(n- 1)! 


1 ’ 


“ A “ A 


.r\^) I A 

^ {n - 1)! ” ^ 


z = 2,..., n—1, 


( 2 . 20 ) 

where A = det P, and A', A" (j = 1,..., n—1) are the determinants of 
the matrix P, in which instead of its the j-th the columns, vectors d', d" 
are inserted. By virtue of condition rank(P) = rank (P, yo) = n—\ and 
relations (12.201) describe all the solutions of the system (12.191) . 

Therefore, the next lemma follows. 

Lemma 2.4. Matrix 


A (n-1)! ) 

2^ 


c= 


n(n+l) 


1 

2- 3 
1 

3- 4 


{n+l)(n+2) 


n{n+l) 

1 

(n+l)(n+2) 

1 

(2n-l)2n 


( 2 n—l)2nC2n-2 

( 2 . 21 ) 
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and the vector 


a = \ - 


n{n + 1) 


1 

"A 




2 ] •> • ■ • ’ 


..., (-l)-2(„ _ 2)!l + A^,) , 

give all the solutions of equation i2.16\) : moreover for 


C2n-2 >0, — Aj 




n—1 


(n- 1)! 


> Co, 


where Co is a root of the equation 


Co 

1 

2-3 


1 

2- 3 
1 

3- 4 


ra(n+l) (n+l)(n+2) 

matrix C is a positive definite one. 

Proof. Let C 2 n -2 > 0, and consider matrix 


1 

n(n+l) 

(n+l)(Ti+2) 

1 

(2n-l)2n 


= 0 , 


Ciz) = 


( ^ 

1 

2-3 


1 

2- 3 
1 

3- 4 


n(n+l) 

(n+l)(n+2) 


\ 


V 


( 2 . 22 ) 


(2.23) 


(2.24) 


(2n — l) 2 nc 2 n- 2 - 


/ 


n(n+l) (n+l)(n+2) (2n—l)2n 

This matrix is positive definite for z = 1/2, by virtue of ii) of the Lemma 
[Q its principal minors A^i, s = arej)Ositive: the Sylvester 

criterion holds. Because of Lemma 12.21 matrix C{z) is positive definite 
for z > fo, where Co is a root of the equation det C(Co) =0. From where 
we obtain that for the parameters C 2 n -2 Cn, which satisfy condition 
(12.231) . matrix C is positive definite. 


Remark 2.1. Observe that equality I2.24\l is reduced to the form 


(l + n) Co + ^ - 


2n 


n 


1 

n+1 


1 

n+1 


= 0 . 


(2.25) 


For establishing this fact, to each i-th row we add all the remaining rows; 
from the obtained rows and columns of the new determinant we then 
extract common factors. For the cases n = 2, 3,4, 5, 6 and 7 the roots of 
Co is equal to 1/3, 5/12, 9/20, 7/15, 10/21, and 27/56 respectively. 
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Lemma 2.5. Matrix = C — CH — HC for the form 
f ^ (A' , A1 i \ 

A (n-l)! J 2 ■ • ■ n 


= 


V 


1 

n+1 


1 

n+1 


is positive definite for 
1 


C2n-2 > 0, 


A 


(n-l)! 


> 


/ 


1 1 
2^^ 


(2n - l) 2 nc 2„_2 
(2.26) 


Co 


1 


4n 


where .Jo A a root of equation \2.25\) . 


Proof. By i) of Lemma [2.31 for s = n and fc = 1, matrix is positive 
definite if the parameter an satisfies the condition 


“ f A' ^ A" 1 - 1 

A ( A, <^n + I — 1. 


1 

n+1 


(2n — l) 2 nc 2 n -2 


(2.27) 


2 

A (^-1)! 

From Lemma 12.21 for 02^-2 > 0, matrix 

( f ? i ^ 

C\z) = 2 3 

, i' 111 

\ n n+1 

is positive definite for all 

^ > J, (2.28) 

where J is a root of equality det C'^(J) =0. By equality (12.251) . we have 

1 1 


b / 


2n-l 


0 = det Ci(J) = det CM 1 + - Jo + 7 : - 


2n 


and consequently 


c= 1 


Co 


2n 


(2.29) 


Thus, the proof readily follows by comparing (12.261) and (12.271) . as well 
as (|2.29l) and the condition p.28l) . 
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3 Synthesis of the bounded controls for the 
canonical system 

The solution of the synthesis problem for the canonical system (ED when 
the controllability function is time of motion, gives the following solution 

Theorem 3.1. Let 

(3.1) 

where ® '^oot of equation iSi.25\) . Let the matrix C and the vector 

a be as in i2.21\) and i2.22\) . respectively. The number ag satisfies the 
condition 

0 < Oq < , -r; F ^ = DnCDn, (3.2) 

2(L ^a,a) 

and the controllability function 0(a;) for x p 0 is defined by equality 
f2.,9l) . where F = D~^C~^D~^, and 0(0) = 0. 

Thus, control u{x) f2.4D transfers the arbitrary initial point x G R" 
to the origin along the trajectory x = Aqx + bou(x) in time T{x) = 0(x) 
and satisfies the restriction |u(a;)| < d. 

Proof. By Lemma 12.41 and condition (13.IL matrix F~^ = DnCDn, and 
consequently matrix F are positive definite. 

If matrix [F—FFI—FIF) is positive definite, equation (12.31) tor x p 0 
has a unique continuously differentiable solution 0(x); see [3]. 

We establish for which parameters a„ and C 2 n -2 is matrix (F—FH—HF) 
positive definite. The positive definiteness of this matrix will follow from 
the positive definiteness of matrix {F~^—HF~^—F~^H). Since matrix 
F~^ = DnCDn, and matrices and FI commute with each other, then 
the following equality is valid 

F-^-HF-^-F-^H = Dn{C - HC - CH)Dn. 

By Lemma 12.51 and condition ()3.1|) . matrix {C—HC—CH) is positive 
definite; hence matrix (F—FH—HF) is also positive definite. 

Therefore, under conditions ED equation (12.31) for x p 0 has a unique 
continuously differentiable solution 0(a;). Set 0(0) = 0, and 0(a;) be¬ 
comes a continuous function for all x. 
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Now let us establish that the control is bounded. We estimate the 
expression a*y{Q, a;)0“ 2 . To this end, we fix 0 and solve the problem of 
finding the extremal of the function a*y{Q,x)Q~^ under the condition 

{Fy{Q, x), 2 /( 0 , x)) - 2ao0 = 0. (3.3) 

This problem we solve with the help of the Lagrange multipliers. The 
Lagrange function has the form 

a* 2/(0, x)0“5 _ X[Fy{Q, x), 2/(0, x)) + 2Aao0. 

Let 2/0 be the extremal point. The necessary condition gives a 0“2 — 
2XFyo = 0, from where we have 2/0 = l/(2A)0“2F“^a. Setting 2/0 in 
restriction ()3.3L we have 1/(2A) = ±y/2ao/{F~'^a,a) 0. Consequently, 
we obtain (a,2/o)0~^ = ±-\/2ao(A“^a, a). Hence, by using the control 
it(x), we obtain 

|w(x)| < \/2ao{F-^a,a). (3.4) 

Selecting the number ap from condition (13.21) . from inequalities (13.41) it 
follows that u{x) satisfies the given restriction in all the phase space. 

Hence from Theorem 1 [3], control u{x) solves the synthesis problem 
of bounded controls in all the phase space, and the time of motion from 
point X to the origin is equal to the function of controllability at x: 
T(x) = 0(x). 


Example. Let us consider the synthesis problem for the system 

Xi = U, X2=Xi, X3=X2, (3.5) 


with a restriction on the control of the form |m| < 1 . 

Since A = —1/20, A( = 1/3600, and A'/ = —1/60, then vector a 
(I2.22L matrix C (I2.2ip and matrix ZI 3 take the form 



Thus, the matrix F ^ = D 3 CD 3 and its inverse matrix has the form 

/ 10- f| -5 I \ / 12 60 120 

F-^= -5 i -I C4, F=--^ 60 180-403 240-12a3 

VI -11/ \ 120 240-1203 -4003 
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From (12.241) and we obtain that = 5/12. Condition (13. 1|) then has 
the form 

„ 1 03 4 

^ ’ 3 ~ ^ ^ 9’ 

according to it, we set C 4 = 1 , 03 = —45. As a result, 


-6 \ 


' 55 

-20 

5 \ 

/ 1/5 

1 

2 

-25 

-20 

10 

-3 

, F = 4 1 

6 

13 

-45 / 

^ \ 

5 

-3 

1 / 

V 2 

13 

30 


We select the number uq of equality (12.31) from conditions (13.21) , which in 
this case has the form 0 < oq < 2/205. We determine the controllability 
function 0(a:) for a: 7 ^ 0 as a positive solution of equation (12.31) (which is 
unique). This equation in this case has the form 


0® = 4104a;?+410030:1X2+82002x1X3+123002x1+53300X2X3+61500:^. 

(3.6) 

The control u(x) ( 1 ^ . which solves the global synthesis problem for sys¬ 
tem (13.51) and satisfies the restriction |u(x)| < 1, is given by the formula 


. ^ 6 25 45 

“<*> = 

Let us hnd the trajectory of the system (|3.5I) . which corresponds to 
the control u = u{x) of form (liTfl) and starts at the point x( 0 ) = Xq G 
This trajectory represents a solution of the system 


25 


45 


Xi = 


0 (x)^^ 0 ^(x)'"^ 03(x) 


X2 


7 X 3 , 


X 2 = Xi, 


X 3 = X 2 . 


Since 0(x) is the time motion from point x to the origin, that is, the 
equality 0 (x(t)) = — 1 is satisfied, thus, 0 (x(t)) = 0 o — t, where 0 o is 
a positive root of equation (13.61) for x = xq . Consequently, the trajectory 
of the solution is the solution of the Cauchy problem of the form 

6 25 45 . . 

xi = - j 7 -72^1-77^-1123^2-777-7732:3, a:2=xi, X 3 =X 2 , 

00 - 1 (00 - ty (00 - ty 

Xi(0)=X?, X2(0)=xl, X3(0)=X^. 


This equation is reduced to the differential equation of the form 

(00 — t)^ Xg^^ + 6 (00 — t)^ X 3 + 25 (00 — t)x3 + 45x3 = 0, 

with initial conditions X 3 ( 0 ) = x®, ± 3 ( 6 ) = x®, and X 3 ( 0 ) = x®. With 
the change of variables t = Qq — e'^, this Euler-type differential equation 
is reduced to a differential equation with constant coefficients with 
respect to 
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the function 2/(t) = X 3 ( 0 o — e"^), which has the form y'" — 9y" + 
33y' — 452 / = 0- Thus, we have 

2 /(r) = ^ci + C 2 cos Vg t + C 3 sin Vg , 

where the constants ci, C 2 , and C 3 are found from the conditions 

y(To) = xl, y'( tq) =-&0X2, y"( tq) - y'( tq) = elxl, (To = ln 0 o) 


and are equal to 


Here 



6 = 


1 


9 0 ^ 5 I 


^0 


6 = - 


VG 00 \0o 



Since a; 3 (<) = 2/(ln(0o — t)), and the functions X 2 it) and a;i(t) are found 
by taking the derivative of the function 13 (t), 


x{t) 


( (00—t) ( 6 ci + G'JG ^2 cos a{t) — by/G sin Q:(t)) 

(00-t)^ (-3ci- ( 3 ^ 1 +v^ ^ 2 ) cosa(t)+ (-3^2+'/6 Ci) sina(t)) 
\ (00 - t)^ (ci + cosa(t) + ^2 sina(t)) 


where a(t) = \/6 in (1 — t/Og) ■ Obviously, x(t) —^ 0 for t —>• Oq. 

The control u(x) on the trajectory x(t) has the form 

u(x(t)) = — 6 ci + 5 ( 6^1 — ^/G ( 2 ) cosa(t) + 5('/6 + 6 .^ 2 ) sina(t). 


For simplicity we solve the problem of attaining the origin from 
the points belonging to the curve Xi > 0 , X 2 = —41a;f/121, and 
X 3 = 0. In this case, from (EH), the time of motion 0 o from Xg = 
(x?, —41(a;?)^/121,0)* is equal to dlx?/!!. The trajectory which begins 
in this point has the form 


x(t) 


(b + 5cosa(t) + 5-\/6 sina(t)^ ^ 
(-6 + 4 cos a(t) - sVG sin a(t)) 

V - 327426 - cos a(t) + Vg sin a(t) j J 
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where a{t) = ^/6 ln(l — llt/{Alx\)). The control u{x) on this trajec¬ 
tory is determined by the relation 

u{t) = (6 + 35cosa(t)). 

Clearly, this control satisfies the restriction. 
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